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ABSTRACT 

The anomalous magnetic moment of the top quark may be measured 
during the first run of the LHC at 7 TeV. For these measurements, it will 
be useful to have available tree amplitudes with it and arbitrarily many 
photons and gluons, including both QED and color anomalous magnetic 
moments. In this paper, we present a method for computing these ampli- 
tudes using the Britto-Cachazo-Feng-Witten (BCFW) recursion formula. 
Because we deal with an effective theory with higher-dimension couplings, 
there are roadblocks to a direct computation with the BCFW method. We 
evade these by using an auxiliary scalar theory to compute a subset of the 
amplitudes. 
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1 Introduction 



The first run of the Large Hadron Collider (LHC) at 7 TeV promises to yield 
a wealth of data and could lead to hints at physics beyond the Standard Model. 
While we do not know what questions the LHC will answer in regards to electroweak 
symmetry breaking, supersymmetry, dark matter or other new physics, we can be 
sure that during the first run our knowledge of the Standard Model particles will 
increase. In particular, the number of top quarks that will be produced at the LHC 
will be comparable to that produced so far at the Fermilab Tevatron, and a much 
greater sample will be produced at large it masses. This will give us an opportunity 
to probe for interactions of the top quark that might indicate its composite structure 
or coupling to new forces. 

One aspect of this study will be the search for anomalous magnetic moment cou- 
plings of the top quark. The consequences of anomalous magnetic moments of the 
top quark have been considered previously, beginning with the work of Atwood, Ka- 
gan, and Rizzo [1] and Haberl, Nachtmann, and Wilch [2]. These authors analyzed 
top quark pair production; they computed the effect of the color anomalous magnetic 
moment on the total cross section and distributions of the top quarks for this process. 
At the high energies available at the LHC, however, one should also consider the effect 
of radiation of additional gluons. It would be useful to have a calculational method 
that could produce arbitrarily complicated tree amplitudes of this type. 

In this paper, we will discuss a straightforward method for computing ti+ng + m / y 
tree amplitudes of arbitrary complexity. In principle, these amplitudes can be com- 
puted from Feynman diagrams. However, the multiple vertices and the complexity 
of gluon interactions make this a challenge. Already at the level of it + 4g processes, 
corresponding to it production with 2 gluons radiated, there are over 100 Feynman 
diagrams. This number increases greater than factorially with the number of gluons. 
A better solution would be to compute the amplitudes recursively, using either the 
Berends-Giele approach |3j or the more recently proposed on-shell recursion formula 
of Britto, Cachazo, Feng, and Witten (BCFW) Some time ago, Schwinn and 
Weinzierl developed a formalism for massive quarks that uses the BCFW method 
and is computationally very effective for QCD tree amplitudes [5]. 

However, the Schwinn- Weinzierl scheme does not generalize directly to include the 
anomalous magnetic moment couplings. The BCFW method requires good behavior 
of amplitudes as some external momenta are taken to infinity. Thus it is nontrivial 
to apply this method to effective Lagrangians that involve higher- dimension interac- 
tions. Indeed, we find that direct application of the BCFW method is stymied by the 
additional momentum factor in the anomalous magnetic moment vertex. 
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Fortunately, there is a way around this difficulty We find that those amplitudes 
that cannot be computed by direct application of BCFW can be computed using 
an auxiliary theory of a scalar particle that carries the spin internally. Combining 
the results, we produce a compact recursive method. This method introduces what 
we consider a promising approach to the application of on-shell recursion to general 
effective Lagrangians with higher- dimensional interactions. 

The outline of this paper is as follows: In Section 2, we will present our notation 
and review some aspects of BCFW computation. In Section 3, we will analyze the 
use of BCFW recursion for fermions with anomalous magnetic moment couplings. 
In Section 4, we will present a useful rewriting of this theory as an auxiliary scalar 
theory. In Section 5, we will present some explicit calculations that check the relation 
of this scalar theory to the original fermion theory. In Section 6, we will present 
our conclusions and compare our approach to other work on the treatment of higher- 
dimension interactions by on-shell methods. 



2 Review of Spinor Helicity and BCFW Recursion 



The goal of this paper will be to present a method for tree-level calculations in 
the theory 



L = * 



Am 



where g is the QCD coupling, a is a fixed constant color anomalous magnetic mo- 
ment, F^ va is the QCD field strength, and = i[y fl , r y u \/2. The same method will 
generalize readily if the theory also includes a QED anomalous magnetic moment 
term 

"2 ea QE D 



SL = 



3 Am 



(2) 



Throughout this paper, we will use the spinor helicity notation, as reviewed ped- 
agogically in [Q\. Instead of using 4- vectors, we will use as fundamental objects the 
spinor products 

(ij) = u_(i)u+(j) [ij]=u + (i)u-(j) . (3) 
associated with lightlike vectors Pi, pj. These objects are antisymmetric and obey 



Kv>i 2 =i[vir = (R+Pi) 2 =^- 

The spinor completeness relation is written in this language as 

P)\p + P]{p =1> ■ 



(4) 



(5) 
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As an example, the polarization vectors of a gauge boson can be written as 



using auxiliary reference spinors r, s. The spinors r, s are arbitrary, corresponding 
to the gauge freedom of the boson. 

The spinor helicity formalism has been extended for use with massive fermions by 
Schwinn and Weinzierl (SW) [5]. For a massless fermion, the helicity states are phys- 
ically distinct and Lorentz-invariant. For massive fermions, there is no unambiguous 
specification of spin state. In the formalism of SW, a lightlike reference vector r is 
used to specify the spin basis to be used. Starting with the massive 4-vector p, one 
defines a lightlike 4-vector p^ by 

Then the u + (p), u_(p) spinors for a massive fermion are 

(p + m)r) [p + m)r] 

u +xP) = — nr^ — u -xP) = — nr^ — • (8) 

It is straightforward to check that these spinors satisfy the required completeness 
relation. 

We will express the values of fermion amplitudes by taking all fermions and an- 
tifermions to be outgoing. With this prescription, outgoing fermions are described by 
spinors u(p) given by 

(r(p + m) [r(p + m) 

Outgoing antifermions are described by spinors v(p) given by 

(p-m)r) _ (p-m)r] , . 

forfe: iAT' ioiqL: ^r- (10) 



To study the effects of top quark polarization, it is useful to be able to compute 
massive fermion amplitudes for an arbitrary choice of the reference vector r for each 
fermion. We will try to retain that freedom in our analysis. 

BCFW [1] proposed a method for computing QCD amplitudes based on the idea of 
deforming the external momenta by a complex parameter z such that total momentum 
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remains conserved and all particles remain on-shell. The explicit deformation they 
proposed chooses two particles i, j and modifies their momenta according to 

Pi-^Pi-zq, Pj-tpj + zq. (11) 

To keep particles % and j on-shell, q must be light-like and satisfy q ■ Pi = q ■ pj = 0. 
For massless i and j this can be expressed as a deformation of the individual spinor 
components: 

i) -> i) - z j) , i\ i] , 

J')^J), (12) 



At tree level, the deformed amplitude A(z) has only simple poles in z from Feyn- 
man propagators going on-shell. BCFW then considered the object 

dz 



/f^)» (13) 



where the contour encircles z = and is taken to oo. If A(z) — > as z — > oo, the 
integral receives no contribution from the contour at oo and the integral vanishes. By 
Cauchy's theorem, this is the sum of residues of poles in the contour. Then, 

A(0) = -^Res^(z) . (14) 



The quantity on the left-hand side of (14) is the original amplitude to be evaluated. 



The residues on the right occur when the deformed momentum that flows through 
a propagator goes on-shell. This relates lower point on-shell amplitudes to the am- 
plitude of interest. BCFW thus obtain a recursion formula that allows the original 
amplitude A to be computed in terms of lower-point amplitudes. 

More explicitly, the recursion relation is 

iA= y £iA L $) j^iA K Q) . (15) 

cuts L 

The sum runs over cuts through a single propagator that divide the amplitude into 
two parts, with the external leg i in the left-hand amplitude Al{i) and the external 
leg j in the right-hand amplitude An(j). These amplitudes are computed with all 
momenta on-shell and with % and j set to their shifted values. The identity requires 
good large z behavior of the amplitude A(z). If this amplitude does not tend to zero 
as 2 — > oo, extra terms appear from the contour at oo that invalidate the simple 
recursion. 
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3 Large z behavior 



Since the BCFW recursion formula depends on good behavior of the shifted am- 
plitude as z — > oo, there is a danger that the BCFW method will not be valid for 
effective theories that contain non-renormalizable operators. In this section, we will 
show that this is a problem for the model ([!]). Specifically, we will show that tree 
amplitudes in the theory ([T]) can be computed in terms of amplitudes with all + or 
— gluon helicities. However, this still leaves a gap that needs to be filled before all 
amplitudes can be computed from simple components. 

In our analysis of the theory 0, we will only consider shifts of gluon momenta. 
In [5], SW give a prescription for shifting the momenta of external massive fermions. 
However, this analysis works only for specific choices of the reference vector r in ([7]), 
while we would like to maintain the freedom to work with an arbitrary choice of r. 



An arbitrary shift on gluons will, according to ( 12 ), involve an external momentum 
i with its angle bracket shifted and an external momentum j with its square bracket 
shifted. There are four possible helicity combinations of the ij to consider: ++, H — , 

and — h For standard QCD with a = 0, the first three shifts give good z — > oo 

behavior while the last case — h does not allow BCFW recursion. Still, for any two 
gluons, there is an allowed shift, and so any qq + ng amplitude can be reduced to 
3-point functions by BCFW recursion. 

Now consider adding to the theory the anomalous magnetic moment vertex. If 
the shift momentum zq flows into the quark line through this vertex, the vertex is 
proportional to z at large z. A fermion propagator carrying the shift momentum 
behaves as z°, and all other fermion vertices — including the magnetic moment vertex 
with an 0(1) external momentum — scale as z°. If the gluon from the magnetic 
moment vertex is connected to external gluon lines through a tree of gluons, each 
propagator in this tree carrying the shift momentum scales as 2" 1 and each vertex is 
at worst z x . Then, finally, the worst possible behavior of amplitudes as z — > oo is z 2 , 
times the ^-dependence of the external gluon polarization vectors. 

If we take q as the reference vector for the polarization vectors of the shifted 
gluons, these polarization vectors scale as 

V2(<?2) z V2(qj) 

£ - (!) = -^m~ z ' ^ = -vssr ; ■ (16) 

We conclude that, in the three cases of shifts allowed in standard QCD, the shifted 
amplitudes behave at worst as 
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i j large z 

+ + ? 

+ - 1 

- - z 2 



However, the true situation is slightly better. For an anomalous magnetic moment 
vertex that stands in front of a fermion propagator carrying the shift momentum, 

™w . , ^ z 4+ ft' + m 

get? (*9» + *») {zq + k y- m " (U) 

we can rewrite 

Wq v = i{<f 4-q»). (18) 

Since q 2 = ($) 2 = 0, the term cancels the leading term in the propagator, and the 
q^ term either vanishes when dotted into a polarization vector or dots with a q in 
a 3-gluon vertex and thus cancels the leading z term in this vertex. For a magnetic 
moment vertex behind a fermion propagator carrying the shift momentum, a similar 
manipulation applies. This reduces the estimates in the table by at least a factor of 
z~ x . In this way, we see that the H — shift allows a BCFW reduction, while the ++ 
and shifts still may not. 

To resolve these last cases, it is simplest to directly compute the amplitudes for 
2 gluons with one magnetic moment vertex in the case of massless fermions. For 
massless fermions, ([T| is not well-defined. However, a massless fermion can have an 
anomalous magnetic moment, and so we replace m in the denominator of the last 
term in Q with some high scale M. This prescription for massless fermions will also 
be used in the discussion in the Appendix. In standard QCD, the massless fermion 
amplitudes with two + or — helicity gluons vanish. With nonzero a, this is no longer 
the case. We find 

This expression behaves as z° after a ++ shift on the gluons. In contrast 

At - - + — — \ 9 2 a (ql} 2 (ql) 2 , . 

A{q ) = -2M ( 9 l)(12)(2g) • (20) 



This behaves as z~ x after a H — shift. This confirms that our current estimates are, 
in general, the best possible. The BCFW recursion can be used to reduce amplitudes 
for which a H — shift is possible, but, for a^O, it cannot be used in the cases of ++ 
and shifts. 

Using H — shifts only, we can reduce any amplitude for qq + ng to amplitudes that 
involve all + helicity gluons or all — helicity gluons. However, we cannot, in general, 
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go further. We need another method to compute these cases, which are required as 
input to the general qq + ng amplitude. 



For the case of gluons coupling to massless quarks, we have obtained an explicit 
formula for the amplitudes with all + helicity gluons. This is presented in Appendix 
A. We have not succeeded in generalizing this to the case of massive fermions relevant 
for top quark physics. In the next section, we will take up another approach to this 
problem. 



4 An Auxiliary Scalar Theory 

We can make progress toward the computation of the all + helicity gluon ampli- 
tudes by breaking up ([TJ into chiral components and rearranging it into a second-order 
Lagrangian. 

Let (ip,ip) be the left- and right-handed spinor components of so that 




(21) 



In this basis, the Dirac matrices take the form 




(22) 



where cr M 



(1, &y and a" = (1, -By, and 




(23) 



where 




The Lagrangian (II]) becomes 




Now formally integrate out i/j and i\r- . This gives 



i — 
2m 



1 



(26) 



L = ip —m + % 



m - i(ga/2m)a flu Ft lua t 



a 



W ■ D 7p . 
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After Taylor expanding the denominator and using the properties of the sigma ma- 
trices, this becomes 



m 



-D 2 -m 2 -i 



99l 



n=l 



i> . (27) 



In this equation, the factor a in front of a- F in (26) has combined with a term arising 



from the commutator of covariant derivatives to produce the factor 

g L = 2 + a 



(2f 



Thus, we obtain a second-order equation with a term close to the full magnetic mo- 



ment of the fermion appearing explicitly. Note that (28) differs from the Lande g 



factor of the fermion, which would be (2 + 2a). The Lande g factor refers to the 
behavior of the fermion in a background magnetic field. The a ■ B term gets contri- 
butions both from a ■ F and a ■ F . Thus, the n = 1 term in the sum in the last term 
also contributes to the Lande g factor, supplying the missing contribution of a. 

If we had chosen instead to integrate out ip and ip , we would have obtained 
the same second-order action with the positions of a^ v and a^ u interchanged. The 
significance of this exchange will be made clear below. 

The procedure of integrating out components of the quark field is used in other 
contexts in infinite momentum frame quantization [7], light cone QCD [8], and soft 
and collinear effective field theory [9]. For the application here, we would like to 



emphasize that this integration out introduces no approximations. From (27), we 
are able to reconstruct any amplitude in the original theory. Although our new 
Lagrangian is not the most convenient way to obtain the scattering amplitudes in the 
limit m — > 0, it does give the correct answers in this limit, as we will illustrate in 
Section 5. 



To analyze the consequences of (27), it is tempting to drop the series of terms 



with a ■ F and approximate this theory by 



V 

m 



-D 2 - m 2 + i^a^F^ a t a 



(29) 



This theory resembles a relativistic theory of a scalar field, except that this scalar 
retains a 2-component internal spin variable on which a^ u acts. In the following, we 
will refer to this model as a scalar theory even though it does describe spin |. To 
better understand the relation of this theory to the original Dirac theory, note that 
if we start from the Dirac equation with a = 



{ip-m)V = 



(30) 
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and multiply by (i p + m) on the left, we obtain 

{-D 2 - m 2 + |s^F^ a t a )^ = 



(31) 



in which the top two components are precisely the equation of motion from (29) 
with gi = 2. The equation (31) arises in calculating of the determinant of the Dirac 



operator, for example, in the background-field derivation of the QCD beta function. 



In general, there is no justification for approximating (27) by (29). However, we 



are interested here in computing the amplitude for qq plus gluons with + helicity 
only. A configuration of gluons with all + helicities is a self-dual Yang-Mills field [TD] . 
The operator a flu F fJ ' u projects onto self-dual field configurations. Conversely, ~& ilv F ilv 
projects onto anti-self-dual configurations and is zero in a self-dual background [TT] . 
So, precisely for the situation of computing an amplitude with all + helicity gluons, 
we may use (29) as a replacement for (27), with is equivalent to Q. The same 
argument implies that, for computing amplitudes with all — helicity gluons, we may 
use the second-order Lagrangian obtained by integrating out if) and if), which has the 
form of (29) with replaced by a^ u . 



The Feynman rules for the theory (29) are the same as those for scalar QCD, 



augmented with the new vertices from the magnetic moment term. These vertices 
contain 2x2 sigma matrices which must be evaluated in the correct external spin 
states. To compute Feynman diagrams in this theory, we first evaluate the diagrams 
as in a scalar theory with an internal spin. The sum of diagrams will contain a 
product of a^ v matrices. We must then take the matrix element of this product using 
the 2-component spinor corresponding to the components of u or v in @, (JlO]) that 
have not been integrated out. Specifically, to compute the amplitude for an outgoing 
fermion with momentum p, in a spin basis described by the reference vector r, we use 
the spinors 

for q R : [p b , for q L : j^j , (32) 
where p b is defined by a7«. Similarly, for an outgoing antifermion, we use 



for q R : p b 



for q L 



—mr] 
[p b r] 



(33) 



A separate reference vector can be used for each external momentum. Finally, to 



account the factor (1/m) in front of (29), the entire amplitude should be multiplied 
by (1/m). 

At any point, we can break up the products of o^ v matrices using the following 
completeness relation: Let a, h be any lightlike vectors that are not collinear. Then 

a] [b — b] [a 



[ba] 



(34) 
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n+2 



Figure 1: Notation for the qq + ng color-ordered amplitude. 



The object on the right is the identity when acting on a] and b], which are independent 
2-component vectors, so it must be the identity in general. This formula is useful to 
write the right-hand side of the BCFW identity as a pair of amplitudes in the scalar 



theory. On a cut line, the scalar 4-momentum will be on-shell, but (32), (33) will be 
replaced by the more arbitrary spinors a] and b]. 

For this construction to be truly useful, we need to show that it allows a BCFW 
recursion that computes the scalar amplitudes in the case of all + helicity gluons. 
For this scalar theory, we can show that the H — gluon shift is always allowed by 
naive power counting alone. First, if the shifted gluons are only separated by gluon 
propagators, the good large z behavior is guaranteed by arguments from pure QCD. 
The only issues arise if the deformed momentum flows through a scalar propagator. 
To argue that even in these cases, the large z behavior is good, we note that scalar 
propagators scale as 1/z and that trees of gluons which contain one of the shifted 
gluons also scale as 1/z. The usual scalar QCD vertices scale at worst as z for large 
z. Thus, the H — gluon shift is always allowed in scalar QCD. The new, helicity 
violating vertices also scale at worst as z for large z, since they are also proportional 
to momentum, and so the H — gluon shift is always allowed in this theory. However, 
this argument fails for the ++ gluon shift and so the scalar theory has the same 
problem that we found in the fermion theory. 

However, there is another possible shift in the scalar theory that is allowed in the 
case of all + gluon amplitudes. This involves a gluon and an external scalar. Consider 
the qq + ng amplitude with the external legs numbered as in Fig. [1} When the gluon 
3 is shifted by 

3) -> 3) - z h) , (35) 
where l\ is lightlike, and the external scalar 1 is shifted by 

l-»l + z3](Zi, (36) 

we will show that the all + amplitude scales as 1/z for large z. Particle 1 is massive, 
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Figure 2: Two classes of diagrams needed to analyze the large z behavior of the 13 shift of 
the amplitude shown in Fig. [TJ 



and l\ must be denned so that the shifted vector 1 remains on mass shell. This 
requires 

i 2 = (1 + z 3] (k) 2 = l 2 + z[3 1 h) = l 2 ; (37) 

that is, 

[31/i)=0. (38) 
A general solution to this condition is 



m 



= 1 - |3l3> 3 ' (39 > 

that is, that li is the flatted vector of 1 computed using 3 as a reference vector. 

In the fermion theory, this shift could also be defined, but it would restrict the 
choice of the reference vector for the external fermion to r = 3 jS]. In the scalar 
theory, we are free to make this choice of a shift without any restriction on the 



reference vector that will appear in (32). 



We will now prove that the above shift on the legs 1 and 3 has the good z — > oo 
behavior that we have claimed. To do this, we consider the possible forms of Feynman 
diagram that can contribute to an amplitude in this theory with arbitrary numbers of 
helicity violating vertices. We need only consider the left-most piece of the diagram 
that contains the scalar line 1 and the gluons 2 and 3. For this, there are two 
possibilities: either gluon 2 is connected directly to the scalar line and gluon 3 is part 
of a tree of gluons with external legs 3 ■ ■ • k, or gluons 2 and 3 are part of the same 
tree of gluons with external legs 2 • • • k. These cases are shown in Figs. |2] (a) and (b). 



To analyze the trees, we need the explicit expresson for the tree in the case in 
which all gluons have + helicity. This expression, for the case in which all gluons have 
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the same reference vector r, is worked out in the Appendix. The result, for external 
gluons j ■ ■ ■ k, is 



aU ' ' J (rj)(j(j + l))---((k-l)k)(kr) ■ l4Uj 

The dangerous term here will be the one that involves the shifted momentum 3 in 
both terms of the numerator. Note that 3 also appears twice in the denominator, so 
this term is only of order z°. 

Now consider the first case, shown in Fig. (2^a). The scalar propagator scales as 
z~ x . If we take the gluon 2 to couple via a magnetic moment vertex, the vertex will be 
of order 1 and the diagram will vanish as z — > oo. Thus, the only dangerous diagram 
is that in which 2 couples by the ordinary scalar vertex. The leading term in this 
diagram is 

, (r!2] i -i3r)(r3 
(r2) (212] (r3){U)---{kr) "' ' 1 ' 

Notice that the numerator of the last term is a matrix. This is a 2 x 2 sigma matrix 



for which we must eventually take the matrix element between the spinors (32) and 



(33). This term scales as z°. 

In the second case, the value of the first part of the diagram is just that of the 



tree ( 40 ) . The dangerous term is 

3r)(r3 



(r2)(23)(34) ••• (Ax) 



(42) 



These two bad pieces contribute to the amplitude at the same order of g and a, 
and so we may add them together. Then an amazing thing happens. The sum is 

:!r)(, ' :! [(2i2](r3)-(ri2](23)]... . (43) 



(r2)(212](23)(r3)(34)---(£x) 
The quantity in brackets is 

[212) (r3) - [21r)(23) = (r2)[213) , (44) 

by the Schouten identity. Since the leading z term in 1 and 3 is proportional to the 
same lightlike vector, this term cancels in the last product. Then the sum of diagrams 
scales as z~ x and the sum has good behavior as z — > oo. This proves our claim that 
the shift on 1 and 3 generates a BCFW recursion formula for the amplitude with all 
+ helicity gluons. 
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We now have an algorithm for computing any qq + ng amplitude for nonzero a. If 
the amplitude contains both + and — helicity gluons, we can apply H — shifts of the 
gluons to reduce the amplitude to lower-point components. If the amplitude has only 
+ helicity gluons, we can use the 13 shift above in the scalar theory to reduce the 
amplitude to lower-point components. If the amplitude has only — helicity gluons, 
we can use the scalar theory with a ■ F. In this theory, a 13 shift that shifts the 
square bracket of 3 reduces the amplitude to lower-point components. Eventually, 
the recursion gives the original amplitude in terms of on-shell three-point amplitudes. 
Though we have given the argument explicitly only for ([l]), the same strategy works 
when the QED anomalous magnetic moment interaction (|2j) is added to the theory. 



5 Calculations in the Scalar Theory 



Although we have shown that the scalar theory described by (27) or (29) can be 



effective for computing qq + ng amplitudes, some aspects of this theory still appear 
odd. Of these, the oddest feature is the factor of 1/m in front of the Lagrangian. 
Some diagrams will then contain factors of 1/m, and one might worry that these 
would generate bad behavior in the limit m — > 0. In this section, we will display 



some amplitudes in the theory (29) that might provide sanity checks on the use of 
that expression. 

First, consider three-point amplitudes. The scalar amplitude as a 2 x 2 matrix is 

(rl2](s32] 



A(l,2 + ,3) 



9_ 
m 



+ —21 [2 

(s31r) +m 2 (sr) 2 J[ . 



(45) 



where r and s are the reference spinors for particles 1 and 3, respectively. Setting 
gi = 2, the fermion amplitudes can be computed by taking matrix elements in (32) 

(sr){r 12] (s32] 

) ' 

(rl2](s32] 



and (|33]). We find 
.4(1 

A(l 
A(l 
.4(1 



2 + ,3" 



gm 



(rl b )(3 b s) (sSlr) +m 2 (sr) 



+ 9-^ 



2 + ,3 



-g 



;i>2] 5 



gm 



(rs)[Vs] 



i- o+ 



2 + ,3" 



[l b 3 b 
[23f 



2 + ,3~ 



(rl b )(s3s] [l b 3 b ] (s31r) +m 2 (sr) 
gm 2 (rs)[r3 b ] (rl2](s32] 
¥W\ ^ (rlr](3 b s) [l b 3 b ] (s31r) + m 2 (sr) ' 
r2][2s] [rs] (rl2](s32] 



-gm 



rl b ] [3 b s] 



gm 



rl b ][3 b s] (s31r) +m 2 (sr) 



(46) 
(47) 
(48) 
(49) 



These expressions are in agreement with explicit QCD calculations. Taking the limit 
m — > 0, these expressions reduce to the familiar three point maximal helicity violating 
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(MHV) amplitudes 

A(l\2\3-) = -g [ ^ , (50) 

A(1-XX) = 9 [ ^ , (51) 

with A{1+, 2+ 3+) = .4(1-, 2~, 3") = 0. 

At four points, there exist two helicity configurations of the gluons that cannot be 
related by parity. These amplitudes can be computed in the scalar theory; we find 

.4.(1,2+, 3+ 4) = / m7 M_ + £h(<£W 12>[32]|2 412][3 
sK ' ' ' ; y (2123) 2m V (2123) (2123) (2123) 2 (212] / 

v ' ' ' ; m (212] (232] 2m (212] [23] J[ V ; 

To compare to fermion amplitudes, we need to take matrix elements of these 2x2 
matrices. For brevity, we will only consider the c/l = 2, m = case. For the case 
with both gluons with + helicity, the massless fermion amplitude with any helicity 

configuration for the fermions must vanish. For massless fermions, the fermion 

projection explicitly vanishes. Multiplying this scalar amplitude by [a on the left and 
b] on the right and simplifying yields 

This indeed vanishes if either 1 or 4 have + helicity. 

In the second case, in which the gluons have opposite helicity, the projection 
should yield the familiar MHV amplitudes at four points. If we choose both particles 
1 and 4 to have + helicity, the projection vanishes by momentum conservation. If 
instead, particles 1 and 4 have opposite helicity, the projection yields 

^■ 2+ - 3 " 4 '»=^ ( 12)SS41> - <54) 

(55) 

which agree with the standard results. 

As discussed in the previous section, the all + helicity amplitudes are completely 



described by this theory. In fact, from the amplitude in (52), one can verify that, 
order by order in a, this expression agrees with that calculated using from (jlj. How- 
ever, a simple observation on the opposite helicity amplitude in (52) shows that this 
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amplitude cannot reproduce the full result from Q. The result above contains only 
terms proportional to a and a 1 , while the exact answer would also contain a term 
proportional to a 2 . This discrepancy is expected, and it is not troublesome for us, 
since this amplitude in the original theory can be constructed using BCFW directly. 



For g L 
states. 



for 



There is one more interesting cross check that we have made of the form of (|29|) 
= 2, a — 0, and so (29) gives an exact description of (1 
At the same time 



for a = 0, all amplitudes of (111) can 
BCFW shifts on gluons with the helicity combinations H — , ++, — • 
compute every qq + ng amplitude in two ways, first, from ([!]) using 
and, second, from (29), using the 13 shift described in the previous 



checked equality numerically to 6 significant figures for all of these 
n = 8 gluons. 



all gluon helicity 
be computed by 
-. Thus, one can 
gluon shifts only 
section. We have 
amplitudes up to 



6 Conclusion 

We have shown that the BCFW recursion relations can be used to compute all 
amplitudes in a theory with an anomalous magnetic moment. The prescription for 
using BCFW is as follows: 

1. If an amplitude contains at least one — and one + helicity gluon, use the H — 
shift to compute amplitudes in the theory defined by 

C = V(ip- m)tf + — yXuuF^t^ . (56) 

2m 

2. If an amplitude contains only + helicity gluons, shift on the scalar and a non- 
adjacent gluon to compute amplitudes in the theory defined by 



m 



-D 2 - m 2 + i^cr^F^H 



1> ■ (57) 



.2 _2 , ,99l 
2 

The gluon momentum should be shifted in the angle bracket. To compute the 
amplitude with external fermions, project onto the fermion line by multiplying 
by the appropriate wavefunctions on the left and right. 

If an amplitude contains only — helicity gluons, shift on the scalar and a non- 
adjacent gluon to compute amplitudes in the theory defined by 



m 



-D 2 - m 2 + i 9 ^a^F flua t a 



ifi . (5* 



2 

The gluon momentum should be shifted in the square bracket. To compute the 
amplitude with external fermions, project onto the fermion line by multiplying 
by the appropriate wavefunctions on the left and right. 
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We have shown that this is an efficient algorithm for computation of tree amplitudes. 
We hope to present some phenomenological applications of this method soon. 

Our conclusions include the statement that the BCFW recursion formula cannot 
be used to fully construct amplitudes in the original fermion theory. This apparently 
contradicts a result of [12], although in fact the anomalous magnetic moment cou- 
pling falls outside the hypotheses of that paper [13]. More generally, the validity of 
BCFW recursion must be thought through carefully for effective theories with non- 
renormalizable couplings. However, our analysis indicates that remedies for their bad 
large- momentum behavior can be found in some cases. 

A distinct momentum shift useful for studying generic theories was introduced in 
[H]. Instead of only shifting the momenta of two of the particles in an amplitude, 
the authors consider shifting the momentum of all external particles. Explicitly, for 
an amplitude with all massless particles, the shift can be expressed as 

i)^i)+WizX). (59) 

i is any external particle in the amplitude, X is an arbitrary, massless four-vector and 
the coefficients Wi are chosen to conserve momentum: 

£>ii]=0- ( 6 °) 

i 

The dependence on the parameter z is easily determined by considering the dimension 
and helicity constraints on an amplitude in a generic theory. For the case of the shift 



in (59), the amplitude behaves as 

A ~ z s as z — > oo with 2s = 4 — n — c — H, (61) 

where n is the number of external legs, c is the sum of dimensions of coupling con- 
stants in an amplitude and H is the sum of helicities of external particles. An on-shell 
recursion exists when s < 0; in that case, there are more angle brackets in the denom- 
inator of an amplitude than in the numerator. In the anomalous magnetic moment 
theory, this all-leg shift leads to a recursion relation precisely for those amplitudes for 



which BCFW fails. This is easily seen at the four point level from (19) and (20). In 



(19), the amplitude is constructible with this shift because there is one angle bracket 



in the denominator and none in the numerator, while in (20) there is one more angle 
bracket in the numerator and so this amplitude is not constructible. This all-leg shift 
could be another way to compute amplitudes in a theory with an anomalous magnetic 
moment. Unfortunately, its practical use is limited because of the proliferation of cuts 
that one needs to compute. 

Recently, there has been some interest in the literature in finding classes of theories 
in whose amplitudes are constructible using BCFW [T5fT6] . The hope has been that 
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the validity of the BCFW recursion formula would say something about the behavior 
of the theory at high energy. It is remarkable that amplitudes in Einstein gravity 
and, equivalently, in Af = 8 supergravity, are are constructible using BCFW [TTyTB"] . 
It has been hoped that this property is evidence for special simplicity of the M = 8 
theory. Further speculations on this point ought to take into account, one way or 
the other, our result that QCD with an anomalous magnetic moment is also BCFW 
constructible. We hope that the methods discussed here can be used to study other 
realistic or effective theories, and that those investigations will shed more light on the 
high-momentum behavior of non-renormalizable theories. 



A A11-+ gluon helicity amplitude for massless quarks 

For massless particles, we have found an explicit formula for the qq+ng amplitudes 
with all + helicity gluons. The derivation of this formula makes use of the off-shell 
current formalism of Berends and Giele [3 J . The off-shell current with all + helicities 
is needed for other arguments in this paper, in particular, in the analysis of the large 
z behavior of the scalar theory at the end of Section 4. 

We consider an amplitude with a single massless fermion line and n + helicity 
gluons. We would like to compute the correction to the standard QCD result coming 
from the presence of an anomalous magnetic moment. Since the background field with 
only + helicity gluons is self-dual, the a ■ F piece of the magnetic moment operator 
gives zero and only the a ■ F of this operator contributes to the amplitude. This term 
has a matrix element only between a + helicity fermion and a + helicity antifermion. 
All propagators and all other vertices in the diagram are helicity-conserving. This 
means that, for the amplitude to be non-zero, the helicity of both external fermions 
must be + and there must be exactly one insertion of the magnetic moment operator. 

The magnetic moment operator contains both a three-point and a four-point ver- 
tex. Both terms contribute to the Berends-Giele current. The terms simplify, however, 
if we choose all of the + helicity gluons to have the same reference vector r. Our 
analysis here generalizes the results of Berends and Giele [3] obtained for currents 
with standard QCD vertices. 

Consider first the term from the four-point vertex. Computing the first few trees 
emanating from the four-point vertex, we find 



J 4 (l,2) 



= i 




(62) 



J 4 (l,2,3) 



= i 



(1 + 2 + 3) 2 



r)(r . 



(63) 



(rl)(12)(23)(3r) 
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From these expressions, we postulate the general form of this term: 

A " °' = i W2hw »' (64) 

To prove this, note that the BCFW recursion is valid for shifting on any two gluons. 
Thus, to prove the general expression, shift gluons 2 and 3 and use induction. Only 
one term contributes in the BCFW sum and it is given by the above equation. 

Similarly, one can compute the first few all + gluon trees that emanate from the 
three point helicity violating vertex: 

J 3 (l) = -tl][l = -i (rl)(lr) , (65) 
, n - , (l + 2)r)(r(l + 2) + (l + 2) 2 r)(r 

Js(1 ' 2) = ~ % (rl>(12>(2r> ' (66) 

Th o ox , (l + 2 + 3)r)(r(l + 2 + 3) + (l + 2 + 3) 2 r)(r 

J 3 1,2,3 ) = -%— / n/iox/oQX/Q \ • 67 

(rl)(12)(23)(3r) 

These suggest a general form, 

, (! + ■■■ + n)r)(r(l + ■ ■ ■ + n) + (1 + ■ ■ ■ + nf r) (r 

J 3 (l,...,n) = -i , — r . 68 

(rl)(12) ■ ■ ■ \nr) 

and that can again be established by induction. 



Note that the first term in the numerator in (68) gives a matrix element between 
a + helicity fermion and + helicity antifermion. The second term gives a matrix 
element between a — helicity fermion and a — helicity antifermion. However, this 
latter term, proportional to the total mass of the gluons in the tree, cancels neatly 



against (64). Finally, we find 



J a (l, . . . , n) = J 3 (l, . . . , n) + Ml, . . . , n) = -i (1 + ' ' - \ ' ' + W) • (69) 

(rl) (12) ■ ■ ■ (nr) 



This result for the Berends-Giele current of all + helicity gluons is quoted in (40) 
and forms the basis for our analysis at the end of Section 4. 

In the case of massless QCD with a single helicity violating vertex, the ampli- 
tude has the diagrammatic form shown in Fig. [3] We can write the gauge invariant 
amplitude as 

A(q + ,gf,...,g+,q + ) = — J(q;l,...,i) / a {i + 1, . . . , j) /(j + 1, . . . , n; q) . 

0<i<j<n 

(70) 
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Figure 3: Generic form of a qq+ng amplitude for massless fermions with a magnetic moment 
vertex and all + helicity gluons. The cross represents the magnetic moment vertex. 



Here, the factors J are currents with an off-shell fermion leg and any number of 
gluons. These were first found by Berends and Giele in [3]. The factor J a is the 
current from a magnetic moment vertex discussed above. The explicit forms for the 
currents J are 

and 

rr ((j + l) + --- + n + q)r) 

J(j + l,...,n q) = //-iTU r=\ • 72 

(r(j + 1)) • • • {nq) 

The amplitude we wish to compute is the sum over all possible insertions of the 
helicity violating vertex with j — i gluons off of the helicity violating vertex and i and 
n — j to the left and right of the helicity violating vertex, respectively. 

Plugging in the various pieces, the amplitude becomes 

A(a+ a+ a+\ = 9 — V + 1 + ••• + »)((» + 1) + ••• + J» 2 

2M 0< £ j<n (ql)---(ir)(r(i + l))...(jr)(r(j + l))---(nq) ' 

(73) 

where momentum conservation has been used. For n = 2, this equation agrees with 



the expression in (19). This expression is gauge invariant as well. It is important 
to note that gauge invariance follows only after summing over all possible places of 
insertion of the magnetic moment vertex. In massless QCD, this is the end of the 
story. We have explicitly constructed the amplitude with all + helicity gluons and 
on all other amplitudes, one can use BCFW to construct amplitudes. 



Since the amplitude in (73) involves only + helicity gluons, it is also possible 
to look at this amplitude as a solution for the motion of a massless fermion in a 
purely self-dual background field. For such backgrounds, Rosly and Selivanov [19] 
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have developed a special formalism, called the perturbiner method, for amplitude 
computations. In this method, they solve the Yang-Mills equations recursively in the 
number of gluons, then use that solution to evaluate the fermion propagator. Their 
solution can be written as 

for color-ordered gluons, with r the common reference vector for all of the gluons. 
The objects £ n are the solutions to the free equations of motion, 

£ = n t e ikn ' x 

where a n is the nilpotent creation operator and t n is the color matrix. The coefficients 
of the product of £s are the Berends-Giele off-shell currents for all + helicity gluon 



derivation of (69 ). 



configurations. Applying (74) to the magnetic moment vertex gives an alternative 



This analysis becomes much more complex in the case of massive fermions. In the 
massive case, the fermion propagators now have helicity-violating factors, and so we 
can insert any number of magnetic moment vertices into an amplitude. In principle, 
we can still construct the all + gluon amplitude with the stitching procedure used 
in the massless case. However, to do this, we need to know an explicit form for the 



analog of the off-shell current in ( 71 ) for massive fermions. In addition, we would need 
to know this current for both helicities of the massive fermion. We do not show them 
here as their form is not illuminating. However, rather than suggesting a general form 
for this current, the expressions seem to get only more complicated as the number 
of gluons increases. It seems that for massive fermions, this method is not useful for 
determining the amplitude with all + helicity gluons. 
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